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Synchronous imaging is used in dynamic space-domain vibration profile studies of capacitively 
driven, thin nH- doped poly-silicon microbridges oscillating at rf frequencies. Fast and high-resolution 
actuation profile measurements of micromachined resonators are useful when significant device non- 
linearities are present. For example, bridges under compressive stress near the critical Euler value 
often reveal complex dynamics stemming from a state close to the onset of buckling. This leads 
to enhanced sensitivity of the vibration modes to external conditions, such as pressure, temper- 
atures, and chemical composition, the global behavior of which is conveniently evaluated using 
synchronous imaging combined with spectral measurements. We performed an experimental study 
of the effects of high drive amplitude and ambient pressure on the resonant vibration profiles in 
electrically-driven microbridges near critical buckling. Numerical analysis of electrostatically driven 
post-buckled microbridges supports the richness of complex vibration dynamics that are possible in 
such micro-electromechanical devices. 



INTRODUCTION 

Suspended resonant nano- and micro-electro- 
mechanical systems (NEMS and MEMS) find use 
in versatile applications, such as ultra-sensitive mass 
detectors, rf filters, and switching devices As 
device miniaturization advances, optimization of the 
overall characteristics in high-frequency MEMS/NEMS 
resonators becomes increasingly complex and linked 
with various mechanical, electrical, thermal and optical 
parameters of the system and its environment. This 
compounds their seemingly superior sensitivity to envi- 
ronmental conditions, such as the pressure, temperature 
and chemical composition of the surrounding gas. 

In the characterization of NEMS and MEMS under 
periodic electrical actuation, vibration profile (VP) mea- 
surements are important in conjunction with frequency- 
domain spectral studies While the latter yield im- 
portant mechanical properties, the former can be useful 
in many applications, including optimization of the exci- 
tation parameters, aiding the identification of sites most 
effective for localized functionalization to enable sensing, 
and in studies of dissipation effects such as intrinsic and 
pressure-dependent damping [3]. Space-domain profil- 
ing is crucial in the presence of significant nonlinearities 
where boundary conditions become critical [6]. For in- 
stance, in typical capacitive electrical drive configuration, 
the force between the grounded substrate and a device 
fabricated by patterned suspended poly-crystalline sili- 
con (polySi) film (serving as an electrode) is inherently 
nonlinear with the drive amplitude and film stress [l|-|6| . 



With interferometric reflection-mode optical transduc- 
tion [4], thermoacoustic effects can significantly mod- 
ify the effective device stiffness or induce autoparamet- 
ric optical drive We observe all these effects to 

be significant in NEMS/MEMS devices defined on films 
with low compressive residual stress under applied loads 
near to the Euler critical value 0- At the critical 
point, MEMS devices are most sensitive to changes in- 
duced by stress variations in chemically-reactive coatings 
[i, m . Of additional practical interest are possible non- 
uniformities of mechanical and electrical film properties 
across the wafer, originating from growth processes and 
application of anisotropic etch, which directly affect each 
circumferentially-clamped microresonator (/iR) Fast 
space-domain visualization of resonant VPs serves as a 
direct means to study the physics of all these effects on 
the single device micro-scale during its actuation. 

VPs in MEMS are traditionally iniaged optically with 
vibrometric [3] , interferometric [10|, UjJ , or stroboscopic 
I2I . microscopy. Recently, spatiotemporal evalua- 
tion of resonant VPs in high-frequency MEMS /iRs was 
demonstrated using resonant realtime synchronous imag- 
ing (RSI) with a pulsed low duty-cycle nanosecond laser 
[13| . The main feature of stroboscopic imaging is a rapid 
production of time-resolved interference pattern movies 
and static profiles, as well as the fast evaluation of VPs, 
thus supplanting scanned motorized probes that are ex- 
pensive and inherently slow. This technique is applicable 
with mechanical resonant frequencies up to /o ^ IGHz. 
In this paper, we use RSI [13| in averaging mode to 
rapidly characterize the VPs in bridge /iRs close to criti- 
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FIG. 1: (a) Schematic cross section of the devices studied 
and definitions of optical quantities used in the analysis, (b) 
SEM images of bridges of dimensions: 25x6x0.12 fxm^ (left, 
slightly post-buckled), and 20x1x0.14 /im^ (right, flat), both 
with ^220 nm static elevations. 



cal stress as a function of the drive amplitude and ambi- 
ent pressure. The effects of high drive nonlinear ity and 
air damping on the resonator VPs are directly monitored. 



EXPERIMENTAL METHOD 

/iRs are fabricated by standard top-down microma- 
chining methods, where bridges are defined photolitho- 
graphically on compressively-stressed n+ doped polySi 
films, deposited by low pressure chemical vapor deposi- 
tion over a sacrificial oxide layer and wet-etch released. 
Upon release of doubly claniped bridges, residual stress 
is relieved through buckling [3, ^ . The devices are driven 
capacitively with the moving /iR serving as an electrode 
and the sihcon substrate serving as a bottom ground elec- 
trode. The inset in Fig. 1(a) shows the schematics of a 
buckled /iR cross section, as well as definitions we use, 
and Fig. 1(b) shows SEM images of our released bridges. 
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In Fig. 2(a) a schematic of our RSI configuration is il- 
lustrated. A dual channel pulse source feeds the /iR and 
optical imaging pulse source (Ao=661.5 nm) in synchrony. 
The collimated illumination at a glancing angle 40° 
is reflected off the /iR and collected by an objective lens 
followed by a 4/ lens pair, the latter of which is used 
for spatial waveform filtering at the Fourier plane with a 
phase mask [15]. The outgoing light is finally imaged on 
a standard CCD camera. Changes in the reflection with 
respect to the static image of the /iR, due to resonant 
motion, are monitored as a function of the rf source fre- 
quency /o, voltage and phase. The pressure within the 
chamber is set with a vacuum pump and venting tubes, 
and monitored via a Pirani gauge. In order to calibrate 
the physical VPs from measured reflection images, an in- 
terferometric analysis is carried out in the out-of-plane 
direction (shown in Fig. 1(a)), as detailed below. Ap- 
plication of a 50% duty-cycle to the imaging pulses (full 
synchronization with the capacitive drive), high in-phase 
sensitivity to average differential actuation amplitudes is 
attained at the expense of lost temporal resolution. 

For calibration of the physical VPs from measured re- 
flectivity images, a Fabry-Perot interferometer multilayer 
analysis, as a function of the total elevation [16], is per- 
formed using knowledge of the static film elevation profile 
do, thickness t, and the refractive indices of the film (n- 
doped polySi, n=3.916) and substrate ns (single crystal 
Si, n=3.834). The reflectance coefficients are calculated 
from the effective reflectivity matrix, assuming nearly 
normal incidence: 
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where 5j = kjdj is the effective phase of layer j and kj 
is the wave number. Denoting Ro{x,y) the reference im- 
age of the static reflection, the contrast signal measured 
during actuation corresponds to: 
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Under full synchronization of the sampling beam with 
the drive frequency and phase, the observed average am- 
plitude < A> dit transverse position (x, y) is then: 



FIG. 2: (Color online) (a) Schematics of the experimental 
setup, (b) Calibration curves for synchronous imaging assum- 
ing a device with film thickness t=138 nm and static midpoint 
elevation do =220 nm. Left: absolute reflection coefficient R. 
Right: differential reflection AR/Rq. With negative values of 
AR, the intensity contrast in the image is negative. 



{A){x,y) = di{x,y) - do{x,y) 



(3) 



With full sampling synchronization and a periodic 
unipolar square wave excitation, the observed average 
amplitude (A) at transverse position (x^y) is (A) = 
Amax/"^' The glancing angle of illumination (see Fig. 
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FIG. 3: (Color online) Drive amplitude dependence of VPs in the fundamental mode of a critically upward buckled resonator, 
(a) Reference image of the static bridge, (b) Frequency domain spectra with low ac actuation voltage and a constant 5 V dc 
bias, (c) Static height profile of the bridge along the Y direction taken from AFM measurements (0 is defined as the height of 
the trench and known film thickness of 140 nm is subtracted on the bridge). (d),(e) Measured synchronous images at /o=4.1 
MHz with different ac amplitudes; (f)-(i) corresponding VPs integrated along Y [in (f),(g), X-profiles], and along X [in (h),(i), 
y-profiles]. Vertical arrows in (f),(g) indicate diminished actuation signals with high excitation. 



1(a)) introduces an additional scaling factor cos 6. The 
peak VP average amplitude, as measured away from the 
static beam position, is then: 

2 

Amax {x,y) = [di {x,y)- do {x,y)] (4) 

cost/ 

giving rise to a normalization factor of 2.83 in our im- 
plementation, with = 45°. A glancing angle also in- 
troduces shadow effects at the resonator edges, which we 
can easily eliminate with appropriate phase masks at the 
Fourier plane (see Fig. 1(a)). We consider only reflection 
variations at the positions of the /iR itself to constitute 
its real VPs. An example of the total and differential 
reflectance curves as function of the total elevations is 
shown in Fig. 2(b). Using Eqs. (l)-(4), together with 
the calibration curve, the average maximum amplitude 
profiles are estimated. Here we will concentrate on char- 
acterizations of the fundamental (lowest) resonant mode. 

RESULTS AND DISCUSSION 

In Fig. 3 we study a /iR of dimensions (25x6x0.14 
/im^) and a midpoint elevation of 220 nm under low 
pressure settings (P<1 Torr). The undriven /iR is al- 
most flat (see Fig. 3(c)) whereas other slightly longer 
devices exhibit noticeable static upward buckling, sug- 



gesting the existence of a compressive force whose mag- 
nitude is close to critical load. Highly buckled resonators 
have been found as hard to drive electrostatically. Fig- 
ure 3(a) shows a static optical image of the unactuated 
device in its initial reference configuration. In Fig. 3(b), 
the frequency response under low-voltage actuation is 
shown. Even with drive amplitudes as low as 45 mV 
and a dc bias of 5 V, we observe the formation of Duff- 
ing nonlinearity [6] and significant spectral broadening, 
with a sudden frequency detuning between 35 and 45 mV 
drive voltage. An AFM measurement of the static bridge 
height profile, in the transverse {Y) direction, is shown 
in Fig. 3(c), and the profile is uniform in the axial {X) 
direction, indicating a shell-like bridge profile. RSI im- 
ages with intermediate and high ac drive voltages, at a 
frequency corresponding to the maximum resonant am- 
plitude, optimal phase and a dc bias of 5 V, are shown 
in Figs. 3(d) and 3(e), respectively. Following the image 
analysis for amplitude calibration, as detailed in the ex- 
perimental section, and integration along the beam width 
(F), the peak VP amplitude X-profiles are shown in Figs. 
3(f) and 3(g), respectively. The peak VP amplitude Y- 
profiles, integrated and averaged along X, are also shown 
in Figs. 3(h) and 3(i). With intermediate drive ampli- 
tudes the VP shapes are as shown in Fig. 3(f) and 3(h). 
With high drive amplitudes, central regions on the beam 
appear to undergo diminished displacement at the orig- 
inal frequency (Fig. 2(g)). However, detuning of the 
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FIG. 4: (Color online) Pressure-dependent studies of VPs in the fundamental mode of a slightly post-buckled resonator, (a) 
Reference image of the static bridge, (b) Frequency domain spectra (inset) and inverse quality factors (dissipation) as a function 
of chamber pressure, under continuous 315 mV ac drive and a 5 V dc bias, (c) Measured interferometric images at /o=2.4 MHz 
as a function of the pressure, and (d) overlaid F-integrated X-profiles of vibration. 



imaging frequency in these cases to values near multiples 
of the fundamental mechanical frequency and the same 
phase settings show some tiny components of vibration 
at these locations. We interpret this observation as re- 
sulting from either nonlinear electromechanical processes 
inducing transfer of energy to higher harmonics at loca- 
tions of high vibration amplitudes on the /iR, or from op- 
tical nonlinearity due to the measured response crossing 
extreme reflection points. In any case, positions with di- 
minished signal, such as the one indicated by the vertical 
arrows in Figs. 3(f), (g) would clearly not be beneficial to 
employ in phase-locked-loop (PLL) sensing applications, 
using this class of /iRs, at this wavelength. Along the Y- 
profiles, slight localization of the motion at the central 
region of the bridge is also observed with high drives. 

Figure 4 shows studies using a narrow microbridge of 
dimensions (25x1x0.12 jam^) and a midpoint elevation 
of 660 nm (t=120 nm and do=660 nm) under varying 
ambient hydrostatic pressure and constant driving con- 
ditions of 1.2 V ac voltage and 5 V dc bias. This bridge 
is slightly buckled in the upward direction, as observed 
in the static optical reference image of Fig. 4(a). Figure 
4(b) shows the dissipation (inverse quality factor, Q~^) 
of the fundamental resonant mode as a function of pres- 
sure, and corresponding spectra (inset). Different pres- 
sure ranges correspond to well-known dominant dissipa- 
tion mechanisms [j, ll^. The total quality factor Q is 
known to approximately scale according to [l^ : 

l/Q = l/Q^nt + aP (5) 

Where Qint is the intrinsic (material) quality factor, a is 
the coefficient of viscous damping and P is the pressure. 



In the data corresponding to Fig. 4(b), a linear fit yields 
Qint = 154 and a = 1.83 x 10"^ [Torr~^] in this /iR. 

In the current experiment we have succeeded in record- 
ing RSI images of sufficient contrast only at pressures 
below the viscous (gas-dominated) regime, namely, cor- 
responding to the intrinsic and molecular regimes in Fig. 
4(b). It is estimated that the most significant limiting 
factors are the low spectral signal-to-noise bandwidth 
(S/N) at low quality factors (below Q ~ 20) combined 
with diminished amplitudes of motion under external air 
damping. Figure 4(c) shows RSI images of the fiR as a 
function of increasing pressure, with a transition from in- 
trinsic to molecular damping. Calibrated maximum am- 
plitude X-profiles, integrated across the beam width (Y"), 
are shown in Fig. 4(d). Increased errors in the VP es- 
timations result from diminished available S/N, giving 
rise to less accurate numerical fits. We consistently find 
that with increasing pressure, the VPs in this fiR be- 
come suppressed around the regions close to the bridge 
overhang (see Fig. 4(d)). This edge suppression effect is 
not observed in repeated experiments under low pressure 
and drive conditions (0.3 V ac voltage and 5 V dc bias), 
that yield available signal-to-noise close to the detection 
limit, with extracted vibration amplitudes comparable 
to the highest pressure case shown here and with more 
pronounced motion near the overhang. 

NUMERICAL MODEL 

The dynamics of a compressively stressed beam are 
described by the equation jl7|, ll8| : 
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FIG. 5: Schematics of the numerical model. 



1.2 V, and the quality factor was (5=7. It is observed 
that the fundamental resonant frequency is 2.45 MHz, 
which is close to the experimentally observed value (Fig. 
4(b)). The corresponding resonant displacement profile, 
averaged over a single period, is shown in Fig. 6(b). 
Small initial imperfection of 0.05 in the initial buckled 
height, corresponding to an excitation of the second anti- 
symmetric mode, was introduced in order to allow non- 
symmetric mode shapes of the beam. Calculations show 
that while the actual beam profiles are dominated by 
the fundamental mode, the resonant deflection profiles 
(i.e., the differences between the initial buckled shapes 




where now E is the Young's modulus of the beam mate- 
rial, I = b X d'^/12 is the moment of inertia of the beam 
cross section, A = b x d is the sectional area, L is the 
beam length, p is the density, and b and d are the thick- 
ness and width of the beam, respectively. In addition, go 
is distance between the ends of the flat side of the beam 
and the electrode (electrostatic gap), eo is the vacuum 
permittivity and V{t) is the time-dependent actuation 
voltage. In accordance with the definitions in Fig. 5, the 
elevation of the beam w{x), as well as the electrostatic 
force f^{x,t) = —eobV'^ /2{go + w)'^^ which is calculated 
using the simplest parallel capacitor approximation for- 
mula, are considered positive upwards. 

Equation (5) has been reduced to the system of cou- 
pled nonlinear ordinary differential equations by means of 
the Galerkin decomposition with linear undamped eigen- 
modes of a straight beam used as base functions. The 
equaations were solved numerically using the ODE45- 
solver implemented in Matlab. The details of the for- 
mulation and numerical approach used for the analysis 
are found in [17] (see also [18]). 

NUMERICAL RESULTS 

Figures 6-8 show numerical solutions of Eq. (5). In 
all cases, the actuation voltage contained both ac and 
dc bias components, such that V{t) = Vdc + VacCos{ujt). 
Zero initial conditions, corresponding to the post-buckled 
configuration of the beam in rest, were used. In all cases. 
Young's modulus E=150 GPa and density p=2300 kg/m^ 
corresponding to polySi were used. 

Figure 6(a) shows the resonant response of the beam 
with dimensions (25x1x0.12 pm^) and electrostatic gap 
^0=660 nm. The axial force was chosen such that the 
midpoint elevation of the beam above its ends was 155 
nm. The driving voltages were Vdc = 5 V and Vac = 
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FIG. 6: (a) Numerical frequency domain resonant response. 
Midpoint deflection of the beam is shown, (b) Deflection pro- 
files (difference between the actual and initial elevations of 
the beam) averaged over a single period. In (b) the operation 
frequency, (jJ=2A5 MHz, is close to the fundamental mode res- 
onant frequency. In both simulations, the dimensions of the 
beam are 25x1x0.12 /xm^; the electrostatic gap is ^o=660 
nm; the initial elevation of the midpoint above the beam's 
ends, due to buckling, is 155 nm; input voltages are Vdc = 5 
V and Vac = 1.2 V, and (a) Q=7; (b) Q=700. 
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and the actual, time dependent shapes of the vibrating 
beam) could be more complex. 

A decrease in the initial separation between the beam 
and the electrode results in an increased contribution of 
higher modes in the resonant VPs. Figure 7 shows time- 
resolved snapshots of the VPs (relative displacements 
from equilibrium) with the same dimensions as in Fig. 
6, but with ^0=220 nm; Figure 8 shows the same vibra- 
tion profile as averaged over a single period. Small initial 
imperfection of 0.05 in the initial buckled height, corre- 
sponding to a contribution of the second buckled mode, 
was again used as an initial condition. Complex displace- 
ment profiles are clearly observed in this case as well. 

CONCLUSION 

Synchronous imaging has been demonstrated as a ro- 
bust method for direct and rapid observations of gradual 
changes in resonant vibration profiles of electromechan- 
ical microresonators under varying conditions of drive 
and ambient pressure. Synchronous imaging can serve 
as a useful tool for studying fundamental processes in 
resonant MEMS/NEMS, as well as for identification of 
favorable device positions most suitable for sensitive 
phase- locked applications, such as sensors, filters and 
switches. Numerical analysis of electrostatically driven 
post-buckled microbridges supports the richness of the 
complex resonant vibrations that are possible in these 
micro-electromechanical systems. 
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FIG. 7: Dynamic snapshots of the deflection proflles (differences between actual and initial elevations of the beam) corresponding 
to different time sections within a single resonant cycle. The beam dimensions are 25x1x0.12 fim.^; the electrostatic gap is 
^0=220 nm; the initial elevation of the midpoint above the beam's ends, due to buckling, is 98 nm; input voltages are Vdc = 1 
V and Vac — 350 mV, and (5=1000. The operation frequency is Ct;=1.495 MHz. Nine (symmetric and skew-symmetric) base 
functions are preserved in the reduced-order model. 
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FIG. 8: Deflection proflles (difference between the actual and 
initial beam elevations) averaged over a single period, corre- 
sponding to the results in Fig. 7 (same parameters). 



